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ABSTRACT 



This paper has a dual purpose. One aim is to study the evolution of coherent states 
in ordinary quantum mechanics. This is done by means of a Hamiltonian approach to 
the evolution of the parameters that define the state. The stability of the solutions is 
studied. The second aim is to apply these techniques to the study of the stability of 
minisuperspace solutions in field theory. For a \(p^ theory we show, both by means of 
perturbation theory and rigorously, by means of theorems of the K.A.M. type, that the 
homogeneous minisuperspace sector is indeed stable for positive values of the parameters 
that define the field theory. 



I. INTRODUCTION 



This paper has a dual purpose. It has its origin in the study of superspaces in field 
theory, that is, the function spaces of solutions of any field theory (even though the name 
superspace originated in the study of the gravitational field). However, we would like to 
emphasize that the techniques we will use here may be applied to the study of coherent 
states in ordinary quantum mechanics, and the examples we will give are, in fact, equivalent 
to one- and two-dimensional nonrelativistic quantum mechanics. 

The simplest example of a superspace is that of a one-dimensional real scalar field 
ip{z, t). If we expand in a real Fourier series (assuming the domain of (p to be confined 
to — L/2 < z < L/2 with the end points identified) 

/ \ / \ f / \ f 2?i7r2\ / 27171 z \ 1 

ip[z, t) = ipQ{t) + 2^ Ypn{t) COS y-J^ J + ^-n{t) Sm ^ J j : (1-1) 

the evolution of (p (independent of the action that generates the evolution) is nothing more 
than a curve in the space of countably infinite dimension defined by the "coordinates" (fn, 
— oo<n<+oo. Of course, classically any nonlinear action for (p gives an extremely 
complicated infinite set of coupled ODE's for the ipn, and this approach to field theory is 
rarely used in direct calculations. 

Nevertheless, such Fourier expansions were the basis for field quantization for many 
years, and are still used in many contexts. It is not difficult to show that for a real Fourier 
series such as (1.1) the quantum evolution of (f is just the quantum mechanics of a particle 
moving in (^^-space under the infiuence of what may be a very complicated potential (even 
in the simplest cases of nonlinear actions for (p). There have been a number of studies 
of nonlinear field theories which have attempted to glean information about the behavior 
of fields such as (p{z, t) by studying model theories where the configuration space of the 
system is reduced by putting all but a finite number of the equal to zero [1, 2, 3, 4]. 
These are called "minisuperspace" field theories. Notice that such theories are equivalent 
to one-particle quantum mechanics in a space of dimension of the surviving (fn- One may 
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ask whether such theories are just models or whether they could be approximate solutions 
to the full field theory. One method for answering this question is to study approximate 
"coherent states" whose center is supposed to move on a track in (^„-space that is centered 
on some classical path which is restricted to the reduced configuration space. 

This problem leads to the study of the simpler quantum mechanical problem of finding 
consistent approximations to the time-dependent evolution of localized solutions to the 
Schrodinger equation in dimensions corresponding to the reduced configuration space. Of 
course, a problem of this sort is completely independent of any field theory, and we feel 
that our results are useful in the general study of coherent states. 

An interesting method for the study of these states is based on the ideas of average 
variational principles. This approach is widely used in the context of nonlinear waves and 
oscillations [5] and more recently in the problem of soliton propagation under the influence 
of various perturbing effects [6] . The emphasis here will be on obtaining approximate and 
rigorous stability results rather than on detailed calculations of the evolution of the solution 
which is the purpose of many works on wave propagation. These stability arguments are 
of much more interest for field theory, since what has been called the "quantum stability" 
of minisuperspace solutions is directly related to the problem of the use of quantum min- 
isuperspace solutions as approximations, but we will find evolution equations for coherent 
state solutions which could be of use in the study of coherent states in ordinary quantum 
mechanics. 

The technique we will use is based on the consideration of time-dependent trial func- 
tions in the Lagrangian of the Schrodinger equation. The Lagrangian is then averaged over 
the space variables to obtain an effective action which involves only the time-dependent 
parameters of the trial function. The Euler equations of this new Lagrangian give us the 
evolution of the parameters. This procedure gives a consistent way to approximate the 
infinitely many degrees of freedom of the wave function by means of a finite number of 
parameters while preserving the Lagrangian structure. On the other hand, pointwise ap- 
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proximations could, in principle (and often in practice), produce in the truncation spurious 
non-conservative terms. Since we are interested in a Hamiltonian formulation of the prob- 
lem, its stability is not determined by damping terms but rather by the nature of of the 
Hamiltonian. Because of this the approximation must be consistent with the underlying 
Hamiltonian structure. We would like to emphasize once more that this idea gives us a way 
to construct approximate coherent states in quantum mechanics and may be very useful 
in what has become a very active field. 

In Ref. [4] a pointwise approximation to the motion of "Gaussian" wave packets in 
minisuperspaces of a A<^^ theory was presented where the action was 

dzdt, (1.2) 

and an topology was assumed for t = const, slices, identifiying the end points z = 
±L/2. The field if was then expanded in the same Fourier series as (1.1). One- and two- 
dimensional minisuperspaces were studied, with special attention to the quantum stability 
of solutions of the one-dimensional sector as imbedded in the two-dimensional sector. The 
one-dimensional superspace had a state function of the form \[/(x, t), where x was equal to 
(fio of (1.1). This function obeyed a non-relativistic Schrodinger equation of the form 

(1.3) 



f\- 

J 2\dt 



1 2 2 , 4 



with e = X/L. Notice that this is nothing more than a one-particle Schrodinger equation 
for a particular anharmonic potential. 

The behavior of localized packets of the form \E' = was then considered, where S 
was taken to be 

S = a(t)x'^ + f3{t)x^ + -f{t)x^ + a{t)x^+ 

^[x- ^(t)]2 + iB{t)x'' + iC{t)x' + i[P{t) + D{t)]x. (1.4) 

To derive the evolution of 5 a pointwise approximation resulting from inserting ^ = 

in (1.3) was used. The parameters a, P, 7, a, B, C, D were assumed to be of order e, and 

a consistent set of coupled linear ODE's for these parameters was obtained. 



The two-dimensional minisuperspace was one where <^o ^ind one of the (pn were taken 
to be non-zero. This will be discussed in more detail in Section IV. 

The equations for the parameters in both the one- and two-dimensional minisuper- 
spaces could be solved recursively, but the theory predicted secular growth of both the 
"center of mass" of the state and the localization parameter (controling "spreading"). In 
this paper we will use the method outlined above to show that the motions are modulated 
and, hence, stable. Since the motion is Hamiltonian, it is not possible to have asymptotic 
stability, and the averaging results are not conclusive for all time. However, it is possible to 
show rigorously, by applying the Kolmogorov, Arnold, Moser (K.A.M.) theorem [7], that 
special motions are indeed stable for all time. 

The paper is organized as follows. In Section II we study the approximate evolution of 
localized states in ordinary quantum mechanics by means of a variational approximation 
where we use (1.3) as an example. In Section III we present perturbation equations for 
the parameters of the localized state of Sec. II and show that their solutions are bounded, 
and then use the K.A.M. theorem to show stability for the full equations. In Section IV 
we study the two-dimensional minisuperspace field theory using a theorem by Melnikov, 
Poesch, and Kuskin [8] to show that the solutions are stable for a set of full measure in 
the parameter space. Section V contains conclusions and suggestions for further research. 
In Appendix A we present the equations for a pointwise approximation for = , 
using an S similar to (1.4). Appendix B contains the full action for the parameters of the 
two-dimensional minisuperspace solutions. 
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II. VARIATIONAL APPROXIMATION 

In this section we will study the behavior of an approximate coherent state in ordinary 
quantum mehanics for a particle moving in the potential V{x) = (//o/2)a;^ + ex"^, where 
£ = X/L gives us the one-dimensional minisuperspace theory. Since L will be something like 
a "radius of the universe" in the minisuperspace case, we will assume that it is very large, 
and thus concern ourselves with potentials where s is small. In order to find the behavior 
of our state, we will appeal to various approximations, but of course there is nothing to 
say that an exact "coherent state" does not exist for this V{x). In fact, depending on the 
definition one takes for a coherent state, one should be able to find such exact solutions. 

For example, if we take a one-dimensional system moving in an arbitrary potential 
V{x), and assume that a solution of the form \E' = e"*^*^^'*) exists, where S = x — w{t)], 
and define u = x, v = x — w{t), S will satisfy the nonlinear partial differential equation 

iw[w {u - v)]— + -—^ + + 



dv 2 dv? dudv 2 dv"^ 



1 fdS\'^ dSdS 1 fdS^ ^ 



The boundary conditions must be taken such that 5 — > oo as a; — > ±oo. The maximum 
value of will occur where (dS/du) ± (dS/dv) = 0, which can be solved for x = F{t). The 
function w can be chosen to make the peak value follow the classical motion of a particle 
moving in V{x). There should be no a priori reason why (2.1) would be impossible to 
solve for these boundary conditions for many potentials V{x). However, for most V{u) an 
analytic solution for S will be difficult to find, and most solutions for coherent states are 
only approximate. 

The pointwise approximation mentioned in the Introduction consisted in using ^ = 
, with <S from (1.4), in (1.3). However, there are a number of other ways to approximate 
such a state. One that is frequently used in the study of solitonic solutions to nonlinear 
differential equations is to insert a simple Gaussian trial function with time-dependent 



parameters into the action for the theory. In our case we can do the same, inserting such 
a function into the usual action for the Schrodinger equation 

"ti /"+00 

dt / 

'to J —oo 

If we take a Gaussian ansatz for ^, ^ = We~'^ , where W = W{t) and 



rti r+oo 1 

L= dt dx[W*'^t - 7:K^x - V{x)^*^]. (2.2) 

Jtn J -OO 2 



S=^[x- g{t)f + iP{t)x + iC{t)x^ + i(j){t), (2.3) 

the g{t) gives the center of the probabihty packet ^f*^', and ji{t) allows for "breathing" or 
"spreading" of the packet. Plugging this expansion into (2.2) and doing the x- integrations, 
we find 



2 + 



dt (2.4) 



In the S given in (1.4) one had to assume that many of the parameters were of order e so 
that a consistent set of equations could be found. Here we are no longer forced to suppose 
that the deviations of these quantities from their harmonic-oscillator values are small, and 
we will see that a much simpler set of equations emerge. 
The Euler equations for (2.4) are 

(5)=0, (2.5.) 

-g-P-2Cg = 0, (2.56) 
6e 

P-filg- Aeg^ - + 2Cg - AC^g - 2PC = 0, (2.5c) 
1^ 



C 3 1 2 3 

+ - - 7 + 7% - = 0' (2-5e) 



2//2 ^2 ^2^ 4 4^2 2//3 
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together with an equation which results from varying W that gives as a quadrature once 
g, P, fx, and C are known. Equations (2.5) can be solved recursively. Using (2.5b) and 
(2.5d) we obtain 

C = A. (2.6) 

From (2.5c) and (2.5e) we find the equation of motion for the center of the packet, 

6e 

g + ixlg + 4eg^ + —g^0. (2.7) 
1^ 



+ - 2/.2 + 2/.2 - + 24e^2 = 0. (2.^ 



Finally, using (2.6) in (2.5e) we find 

2 

If we introduce the change of variables 11 = we find the final system 

g + lilg + 4eg^ + 6sa'^g = 0, (2.9a) 

a + /ino - A - 12£a^ + 24sag'^ = 0, (2.96) 

for the position g and the variance ( "breathing" ) of our approximate coherent state. 

In order to compare these equations for /j, and g with those in Ref. [4], we can see 
that if we write, as above, /j, — no + 2ej, and g = go + 2ea and keep terms to order s, we 
find 

i?o + /^o^o = 0, (2.10a) 
a + lila^-2gl-^, (2.106) 

and 

7 + 4//2^ = 6 + 12^2- (2.11) 

Instead of comparing these equations directly with those of Ref. [4], it is more useful to 
use the parametrization used in this article and apply it to the pointwise approximation 
of Ref [4] . To do this we write 

S = a{t)x^ + P{t)x^ + g{t)]^ + iB{t)x^ + iC{t)x^ + iP{t)x + i(j){t). (2.12) 



(leaving out a trivial real ^-independent term which appears in Ref. [4]). We have to 
assume that cu, /3, and B are small enough that we only keep terms to linear order in them 
or we will have terms in x of order higher that x'^, and the ansatz will be inconsistent. We 
can substitute this S in the Schrodinger equation (1.3) with ^ = and find a series of 
coupled equations for a, /3, /i, g, B, C, and P (we will ignore the equation for (f)). These 
are given in Appendix A. The equations for a, and B are the same as in Ref. [4]. If we 
take fj, = fiQ + 2^7 and ^ = + 2£cr, and assume that a, P, and B are of order e, we find 
a = 0, a = e/Ahq. Equations (A3) and (A4) imply 

(3 + 9nl(3 - 4eno9o = 0, (2.13) 



and the solutions are (taking go — Q cos not) 



/? = exs cos 3//o^ + cos //o^) (2.15) 



Qs £ 

B = ex^smS/jiot — - — sin/xg^ o' (2.16) 

2//0 3//q 



where ^3 is a constant. Now, Eqs. (A6) and (A8) can be solved for C and P respectively. 
Taking the derivative of the equation for C and inserting the result in (A5) gives 

II 2 \ 1^1^^ 

+ ^g^ + + ^Qf^9/3 + B (^-36^ - 33^^^ = 0. (2.17) 

Substituting the result for P in (A7) gives 



9 + ^^l9 + £[^9^ + + ^99 ) + 



3 , ^9^9 , 6/i2 



+/3(30^^ +12---3^ ---gg- — ] + 
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+s(-36^^,-^V)=0. (2.18) 

These two equations can be compared to (2.7) and (2.8). Note that they are the same 
except for differences in the order-£ "driving terms" that modify the second order ODE's 
that determine g and 

In order to compare these equations with (2.10b) and (2.11), we wiU again take ji = 
Ho + 2ej, g — go + 2ea. To order e we again get (2.10a) for go and 

^ + /^^^ + (3^0^ + ^) + (^) (l5/.o^o^ + 12- 
6^0 A ,.fB 



- 18 - go9o = 0, (2.19) 
/^o y \e J 

7 - + (^3fio9l + ^ + 7^) + IS (^^^ iilgo - 18 (^|^ fio9o = 0. (2.20) 

Once more these differ from (2.10b) and (2.11) in the "driving terms", where the term 
proportional to e is shghtly different and there are added terms in B and /?. 

Clearly, the pointwise approximation, aside from being more complicated, does not 
add much to the simpler equations from the variational approximation. In Ref. [4] the 
main problem with the equations for 7 and a (i.e. h and g) is that they generated terms 
with secular growth, which made the predictions of the theory only valid for a limited 
time. Both (2.17) and (2.18) have terms that drive a and 7 at resonance, so to order e 
in perturbation they will also have secular terms. In Ref. [4] it was assumed that better 
perturbation methods would remove such terms. In the next section we will apply such 
methods to the variational approximation as well as provide a rigorous demonstration that 
the solutions for fx and g are indeed bounded. 
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III. ABSENCE OF SECULAR GROWTH 



A. PERTURBATION SOLUTION 

In the previous section we constructed the action and obtained the equations for /j, 
(i.e. a = and g. Prom (2.10b) and (2.11) one can see that the order-£ perturbation 

theory leads to secular growth for both the g and a perturbations. However, there are 
more sophisticated perturbation methods methods which remove these terms and provide 
a uniform approximation to the solution. 

If we look at Eqs. (2.9) and make certain assumptions about the size of a and g, we 
can readily show that g must be bounded. If in the last term on the LHS of (2.9a) we take 
a to be little enough different from l/y^io so that ea^ fti ^/a^o, (2.9a) then becomes (with 
no assumption about the size of g) 

9=-(f^l + -)9-^eg\ (3.1) 

This is the equation of a particle moving in the potential 

and g can be solved for in terms of elliptic functions. The form of V shows that the motion 
of g is always bounded. Here, of course, we are assuming that a is bounded (and small 
enough) so that there are no large excursions of the 3eg'^/nQ term in the potential. 

In order to show that a is bounded we can take a to be (l/yZ/io) + a, where a is small 
compared to l/^/JIo (but not necessarily of order e). As long as ea is small enough to be 
ignored, we need only take in the last term of Eq. (2.9b). Taking go{t) = Q cos not we 
have that the equation for d is, to second order, 

a + A^ld - 6//o^'a2 = -12 ( ^ ^ 3^/° ) e - 12^Q' cos 2/iot- (3.3) 
To remove the constant forcing to order we take d = p — ^(i+Q mo) ^ ^^^^ obtain the 

Mo 

equation 

p + (4/^0 — eK,p,i)p — Kp^ = eF cos 2//0*, (3.4) 
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where k = 6^1^^, Wi = 6 ( ^^f^ ) , F = -12-^. This equation, for « = 0, is the 
equation for an oscillator forced at resonance. The equation will now be solved to show 
that the nonlinear terms prevent the growth of the solution. To obtain the approximate 
explicit solution we will use the method of averaging for the Lagrangian 



to 



^ ^ 1 .2 ^ M±£^!i)p2 _ _ 2/.ot \ dt, 



(3.5) 



and use the trial function 



p = e^/'^[A{st) cos2/io* + B{st) sm2jiot] + sD{et), 



(3.6) 



where the mean value D has to be introduced due to the quadratic nature of the nonlin- 
earity. Averaging over the fast time 2//ot, we obtain 



Jto L ^ 4 2 J 2 



dt. 
(3.6) 



We now assume the so far arbitrary (but small) amplitude F to be F = s^I'^Fq to balance 
the resonant forcing and the nonlinear effect. With this, e^l'^F fti 0{e^^ becomes of the 
same order as the modulated terms (the argument can, of course, be reversed, assuming 
that F C(l) and rescaling p and the time appropriately). 

Variation with respect to D gives the usual algebraic equation for the mean value as 



2D - -{J^ ^B"^) = 0, 



(3.7) 



I.e. 



D = -(A^ + B^). 
4 

Using this result, and integrating by parts, we obtain the final form of L as 



(3.8) 



L = 2£>o 



to 



AB' 



(A' + B') + ^{A' + BY + 
3z/io 



8/10 



-A 



dt. 



(3.9) 
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Clearly the solutions of the Euler equations are the level curves of the Hamiltonian 

H^-^{A^^ B') + ^{A' + B'f + -^A. (3.10) 

It is easy to see that the critical points of the system are the critical points of H. They 
satisfy dH/dA = dH/dB = 0. the local minima are centers and the local maxima are 
saddles. Also for large A'^ + B^ all curves are closed. It is then clear that all the motions 
are bounded. This shows that the resonant forcing is balanced by the the nonlinearity, 
giving a bounded motion for p. This in turn implies that the coherent state does not 
spread indefinitely away from its original spreading (//q)"^^^. To complete the description 
of the motion we now examine the central points. They satisfy 



dH F fiiA K 



2 



0= -^ + ^^ + — A(A2 + s2), (3.11a) 



dA 4:^0 4/io 8;Uo 

dH 



dB 



= B 



^1 \+J^(A^ + B^) 



(3.116) 



,4/io/ 8/io 

The only solution for the second equation is S = 0, since /ii is always greater than zero. 
When B = the equation for A is the cubic 

^2 

F + fiiA+ —A^ = 0. (3.12) 

There is only one critical point (^*, 0) where A* is the only solution of (3.12). Clearly this 
is a minimum of H. This represents a periodic solution of finite amplitude A* with period 
2//0- In this case all level lines of H are concentric curves. They give periodic solutions for 
A and B, which in turn represent quasiperiodic motions of p. 

This perturbation solution shows that small perturbations remain small. Interestingly 
enough some of these perturbation results can be proved rigorously as we will now show. 

B. RIGOROUS RESULTS 



We now establish the boundedness of the motion for not necessarily small initial 
conditions of the amplitude g of these oscillations. We consider again the Eq. (2.9b). Here 
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we will assume that = Acos/jiot, where A is not small. Moreover, assume the initial 
conditions for a to be arbitrary, that is, not necessarily close to l/yT^. In this case we 
will show that for sufficiently small e the solution for a{t) is always bounded. To do this 
we appeal to the classical results of periodic perturbations of Hamiltonian systems with 
one degree of freedom [7, 9]. This pertubation theory is based on the Poincare map. Since 
the equation is periodic with period 27r/2//o this map is defined as follows. Consider the 
map: 

[a(0),d(0)] ^ [a(27r/2//o),a(27r/2//o)], 

where (a, a) is the solution of (2.9b). Clearly this map is area-preserving since the system is 
Hamiltonian. An invariant curve of this map represents a two-dimensional invariant torus 
in the space (a, a, t); {t mod 27r/2//o)- Since two-dimensional tori divide three-dimensional 
space, we have that initial values inside an invariant torus remain there. This bounds the 
solutions. Thus the problem of boundedness of solutions for all time is tranformed into 
the problem of finding invariant curves for the Poincare map of (2.9b). 

The existence of invariant curves is proved using action-angle variables for Eq. (2.9b). 
To study this map it is convenient to first introduce action-angle variables for the unper- 
turbed oscillator. The transformation is not expressed in terms of elementary functions, 
but an explicit representation is not necessary in order to derive the results. 

The generating function for the canonical formulation is given by 



W{a,E)= l\lE-(^e + ^]dC. (3.13) 







where the variable E is expressed in terms of the action (the area enclosed by the orbit) 

Clearly I'{E) > and E can be given as E{I). The generating function becomes W{a,I). 
We now have 

P=a=^, 9=-^ (3.15) 
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in the new canonical variables / and 9. In these new variables the Hamiltonian is inde- 
pendent of 9. 

The perturbed problem in these variables takes the form 

/ = -e-^(^,/,2//ot), (3.16a) 

e = + ^^^(^, '^l^ot) (3.165) 

The Hamiltonian H^^^ is the image of (a^/12)^^ under the action-angle change of variables. 
The explicit form is not needed, just note that H^^^ is a smooth function of the new 
variables. The Poincare map in these variables, if (/q, ^o) is the inital value and 9^^^) 
the final value, takes the form 

/W=/o-£ ——[e(^j,^eo),Iii,Io,9o),2fioi]dC, (3.17a) 

^^'^ = ^0 + j -^m,Io,eo)]+eJ ——[9{^,i,,9o),m,Io,eo),2fioC]dC 

(3.176) 

The integrated terms are smooth functions of Iq, 9q since the solutions depend smoothly 
on the initial values (/q, ^o) and the integration interval is finite. Thus the map takes the 
form 

=/o + £F(^o,/o), (3.18a) 
^(1) ^9o + —u;{Iq) + eG{9o, /o), (3.186) 

where a'(/o) = ^W/^I\I^I^^. 

Since / can be expanded in a power series in e (the function is analytic in e because 
the equation is integrated over a finite interval). The integrated terms of a;(/) contribute 
to the term G to order e. The contribution to the leading order just comes from Iq. 

This is the canonical form of the map for the application of the K.A.M. theorem 
[10]. The map for e = has as invariant curves the circles / = constant, < ^ < 27r. 
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Moreover on each circle the rate of advance is ttuj{I)]/ ijlq which is action dependent. (In 
this simple case the invariant curves of the map are just the level lines of the unperturbed 
Hamiltonian.) In this setting the K.A.M. theorem guarantees that for E sufficiently small 
and suflBciently irrational the invariant curves persist [10]. It is also known that this 
set of "sufficiently irrational" numbers has full measure. Thus, since the function is 
increasing (and tends to infinity as / — > oo) it is clear that there is a set of full measure in 
the variable / for which invariant curves persist. It then follows that the motion of a(t), 
a{t) is bounded for all times, provided that the initial conditions fall within the perturbed 
invariant tori. This gives a rigorous verification of the perturbation results obtained in the 
previous section. 

Up to this point we have made no reference to field theory. The results in this section 
are actually a convenient way to study approximate coherent states in one-dimensional 
quantum mechanics. They can be applied to the motion of coherent states with "breathing" 
or "spreading" such as those studied by Guth and Pi [11]. There are a number of ways 
to extend our results to field theory, some of which will be mentioned briefiy in Sec. 
V. As mentioned in the Introduction, one of the main objectives of this article was to 
use the variational approach to investigate how close a microsuperspace solution is to a 
larger minisuperspace solution over a long period of time. In order to do this, in the next 
section we will apply the methods used above to the same minisuperspace-microsuperspace 
problem used in Ref . [4] . 
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IV. MINISUPERSPACE FIELD THEORY 



In Ref. [4] coherent states peaked around some small sector of a field theory were used 
to investigate whether they could be approximated by a state restricted to that sector and 
then quantized. As discussed in the Introduction, the system was a Xip"^ model in an 5^ 
topology with only two modes "unfrozen". That is, if in (1.1) we put all but (fo and one 
of the (pn equal to zero, we get a minisuperspace field theory by plugging this ansatz into 
the A<^^ action. Putting the <^_„ equal to zero and dropping third and fourth order cross 
terms among the (fin in the action yielded a two-dimensional classical theory that could 
be quantized. Moreover, using the definitions ipo = x and (pn = y it was shown that the 
minisuperspace state function obeyed the Schrodinger equation 

where s = X/L and ttiq = /Iq + (27rn/L)^. In Ref. [4] a pointwise approximation for 
= e"*^ with S similar to (2.12) was used. Here we will use the variational approach with 

^ = W{t)e~'^, where 

2 2 

Li TYl 

S=^{x-gif + —{y-g2f + exy + iPix + iP2y + iCix^ + iCay^ + iMxy + i0(t) , (4.2) 
and /i, m, gi, g2, 0, Pi, P2, Ci, C2, M are all functions of t. If we insert this S into 

rti roo I'oo i 

L= dt dx dy[i***t--*;*,--*;*j,-y(a;,y)***], (4.3) 

Jtn J -00 J -00 ^ ^ 



with 



V{x, y) = ^x' + ex' + ^y' + 6exV, (4-4) 



we find a very large action for /i, m, gi, g2, 0, Pi, P2, Ci, C2, and M which is given 
in Appendix B. Of course, we are only interested in states that can be approximated by 
microsuperspace states with y — 0, i.e. those with g2 = 0. If we now look at the form 
of ^f*^', we see that a ^'*^' = const, surface is an ellipse centered on y = 0, a; = ^1 
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with its semi-major and semi-minor axes given by and m (which axis is associated with 
which parameter depends on their relative size). The parameter 9 gives the angle that 
this ellipse makes with the x-axis. There is some freedom in defining what one means 
by a minisuperspace state "close" to a microsuperspace state (here the y = state). For 
simplicity we would like to take ^ = and g2 = 0. If one varies the action (Bl) with 
respect to g2, P2, 0, and M, one finds a series of equations which, for the initial conditions 
= ^2 = 0, -P2 = -P2 = 0, ^ = 6* = 0, M = M = 0, maintain 6*, ^2, -P2, and M zero for all 
time. For this solution the reduced action 







(4.5) 



gives the correct equations of motion. 



The Euler equations for this action are 



2C2 



+ 



m 



= 0, 



(4.6a) 



m 



2m? 





-gi - Pi - 2Cigi = 0, 



(4.6c) 



A + 2giCi - iilgi - 4egl - 



6egi 



- 2P1C1 = 0, 



(4.6d) 



-ACfg, - -g, = 0, 



(4.6e) 






(4.6^) 
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together with the normahzation condition {W"^ / ^mjj)' = and an equation for the trivial 
phase (f). These equations give 

-g,-^,lg,-4sgf-'-^-'-^ = 0, (4.7a) 

fi m 

_/i ^ 3 //iy _ 2 2 + 2;.^ ^ 24s^2 + 1^ + 1^ ^ q, (4.76) 

-- + I f-)' - + 2ml + 2% f^? + = 0. (4.7c) 
m 2 \m J \ 2/1 J 

In this case we can also appeal to a K.A.M. type of result in order to obtain qualitative 
information about the nature of the solutions. For this we observe that the unperturbed 
system 

gi + l^ai = 0, (4.8a) 
d + ula - ^ = 0, (4.86) 

b + (4.8c) 

can be written in terms of action-angle variables. The equations for a and 6 have the 
variables given in (3.13) and (3.15). The equation for gi has the polar coordinates as 
action-angle variables. Denoting the actions by Ig, la, lb and the angles by 9g, 9a, Ob, the 
perturbed problem takes the form 

'^ = -'-89;- ^^ = ^°+^^' ^'-'"^ 

^« = -^^B— = M^a) + £^1—, (4-96) 

= -^^E- = Mh) + £^1—, (4.9c) 

Oub oh 

Since the frequency jiQ does not depend on the actions, the usual K.A.M. theorem cannot 

be applied to guarantee the persistence of the invariant tori Ig = const., la = const., lb = 
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const., instead a modification due to Melnikov, Poesch and Kuskin [8] can be used. The 
result is as follows. If the unperturbed Hamiltonian system is of the form 

= A, i = 0, (4.10) 

where A is a constant vector, then small Hamiltonian perturbations preserve the invariant 
tori I = const. = At + for most values of the vector A. In this case we choose 

A=[//o,a;a(/°),a;6(/5°)], (4-11) 

and thus the remainder is small provided that we search for tori close to 1°, /°. The result 
guarantees that for most values (a set of full measure) of hq, I^, the unperturbed tori 
persist. 

This shows the existence of quasiperiodic, and thus bounded, motions for most initial 
conditions In this case the existence of invariant tori does not prove stability 

because three dimensional tori do not separate the six-dimensional phase space. Moreover, 
a universal instabilty, known as Arnold diffusion is present. This implies that for sufiiciently 
long times, that is, tAD — 0{e^/^) particular initial values always leave any bounded region 
of phase space. In general this time is very long compared to other time scales in the 
problem. In order to see whether this instability is relevant for our problem, it is necessary 
to rewrite our equations in terms of more conventional units. The expression for the Arnold 
diffusion time mentioned above is actually tAD = (l/'^o)^^/'^, where e is a constant which 
estimates the ratio of a perturbation to a Hamiltonian to the unperturbed Hamiltonian, 
and u>o is the frequency of oscillation associated with the unperturbed Hamiltonian. 

The units of if in (1.2) are Q£, where £ is the unit of length and Q is a "charge" 
associated with (p. The fact that we are working in one space dimension rather than three 
means that has units of force. The constant /iq has units of one over length, and the 
constant A must have units of the reciprocal of energy times length cubed, and e has units 
of one over energy times length to the fourth power. This, along with the fact that d/dt 
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is really d/dct, allows us to determine the units of g and a which are both Qi^l'^. We can 
define the dimensionless variables, u = g/QL^/"^ and h = ajQL^I'^^ and Eqs. (2.8) become 

u + c^/ZqW + 4£o (^1;) u^ + 6£o(-^) = 0, (4.12a) 



where eo is a small dimensionless parameter. By adjusting Q we can put £o equal to one. 

We need only give /xq and Q to find Iad- These quantities depend on the type of field 
that is supposed to model. Since <p is one-dimensional, these quantities are somewhat 
different from those one would expect in the three-dimensional case. For //q, a classical 
quantity that has the proper units is mc^/e^, where we can take e and m the charge and 
the mass of the proton, which gives loq — mc^/e^. A quantity Q associated with the 
charge that has the proper units is Q = m(? je (note that the dimensionality forces us 
to take a reciprocal of e). This implies that c^Ti' jQ'^L^ is {l/aj){\c/LY{c/L)'^, where 
Ac is the Compton wavelength of the proton and a/ is the fine structure constant. This 
quantity is small with respect to c/L for L > Xc, so Hi is proportional to (c/L)'^, while 
Hq is proportional to m^c^/e^, so tAD is roughly tAD = (e^/mc^) exp[(Lmc^/e^)^]. In 
principle, Arnold diffusion could be important when L < jm(? . However, we can take 
a cosmological model which, for simplicity, has the time behavior of a A; = Robertson- 
Walker universe, but the mass of a A; = 1 model, which gives the age of the universe for 
any L to be 



where Mq = 5 x 10^*^g. 

If we now take L = /?(t)(e^/mc^), then the ratio of tAD to the age of the universe at 
any time would be 

tAo 1 ./3V^![^V^'- J_e/3^(3.5xl02«). 
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This quantity has a minimum of order one at (3 of order one and grows for all other values 
of so Iad is always much larger that the age of the universe. 

If (f is supposed to model the gravitational field itself, we can take = c^/ G (again 
the reciprocal of the usual "charge"). Here c^h'^ /Q^L^ is {Lp/ LY{c/ L)"^ , where Lp is the 
Planck length. Since the field is massless, we can take jiQ = 0. Here, then, the {c/LY 
terms are Hq, and the term in c^f? jQ'^L^ is ffi, so e is of order [Lp/LY, so Iad is roughly 
tAD = {L/c) exp[(L/Lp)^], and the ratio of Iad to the age of the universe (the age of the 
universe does not change) is (now taking L = (3{t)Lp) 

which is, again, always much larger than one. 

We thus conclude that the approximation based on coherent states is good in the 
sense that classical trajectories of the center of mass are not qualitatively modified by the 
existence of small numbers of higher modes in field theory. However, when infinitely many 
modes are present it is not possible to draw any rigorous conclusions about the limit. This 
case has to be investigated independently. Even simple finite-mode systems exhibit new 
features in the limit. 
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V. CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 



The twofold purpose of this paper has been addressed by means of the approximation 
techniques given in Section II. We have shown that it is possible to obtain time-dependent 
consistent approximations to the motion and breathing of coherent states using average 
Lagrangians. These Lagrangians have allowed us to explore not only the evolution of these 
states in quantum mechanics, but have allowed us to study the stability of such solutions in 
minisuperspace field theory. The most severe potential instability arose from the resonant 
interaction of the center of mass of the state and the width parameter. We showed, using 
a technique of averaging of ordinary differential equations, how nonlinear effects saturate 
the resonance amplitude and change the solution into a modulated solution. In special 
cases the stability results can be proved rigorously by means of the K.A.M. theorem. This 
formalism suggests a new potential instability when many modes of comparable size are 
present in the coherent state. In fact, a finite but high dimensional Hamiltonian system 
is obtained. In this case Arnold diffusion could have a destabilizing effect over very long 
times. Clearly the same instability will be present in any Hamiltonian truncation of more 
realistic field theoretic models, and its relevance in practice will have to be studied in each 
case. 

Finally, we would like to emphasize the fact that this variational idea could be applied 
to functional Lagrangians. The parameters will now be functions of both position and time. 
The averaged action will have as its Euler equations nonlinear partial differential equations 
which will describe the spread and the position of the coherent state in function space. 
This is currently under investigation for a simple <y?^ functional model. 
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APPENDIX A. POINTWISE APPROXIMATION EQUATIONS 

If we write = with S given by (2.12), then S obeys 



.dS Id^S 1 fdS^^ 

+ 



dt 2dx'^ 

Inserting (2.12) in this equation and assuming that a^, and B'^ are small enough to be 
ignored, we can equate the real and imaginary coeeficients of each power of x on the LHS 
equal to zero. We find 

x'^ : -Aaji + £ = (real) (^1) 

a - 8aC = (Imag.) {A2) 

x^ : Aajig - Spfj, + 6BC -B = (Real) (^3) 

/3 - AaP - 6Cp - SBfj, = (Imag.) (A4) 

2 2 

x^ : 6a-C-^ + 2C^ + SP/ig + 3PP + ^ = (Real) (Ab) 

^ - 3(3P + ZBjig - 2/xC = (Imag.) (A6) 

x: SP- P + iJ,^g + 2PC = (Real) (^7) 

-fig - iig + 3B + 2^iCg - /iP = (Imag.) (^8) 
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APPENDIX B. COMPLETE TWO-DIMENSIONAL ACTION 



The action (4.3) for = W{t)e ^ with S given by (4.2) is, after integration over x 
and y, 



to yjmpi, — 6'^ 



Piigi- 



m/i — d"^ 



+ P2\g2- 



IjQgi - 9-g2 
m/j, — 0'^ 



+ 



+ 



2{mn - 9'^f 



+ 



+C2 gi + 



29figig2-29^gi , fi9^ - 9ym + 2{^g, - 9g2)^9^ 



2m mil — 9"^ 

,2 



+ 



2(mii - 6*2)2 

( I^rngig2 - n9gl _ 9 _ Ii9{mg2 - 9gif 
V mii-9^ 2{mii-9'^) 2(m// - 6»2)2 

_1 2 fj_ rn9'^ - + 2(mg2 - 9gif9'^ 
2^ \2ii^ 2(m// - ^2)2 



+ 



-Ii9 



lJ-mgig2 - li9gl 
my. — 6*2 



9 



2(m//-^2) (m;U-^2) 



H9{mg2 - 9giY 
2 91 



92 - 



li9gi - 9'^g2 
mn — 6*2 



Jlf 2. 1 2i,9g^g2-29^gi u9^-9ym + 2{iig,-9g2)^9^ 
2 2m mn-9'^ 2{mii - 9^ 

^_nr^ f 2,J__ '2rn9gig2 - 29^ gf m9^ - 9^ fi + 2{mg2 - 9g^)^9^ 
- 2 ^^1 [91 + mii-9'' 2{mg2 - ^2)2 

_iM2 f ^ + 2/^%i^2 - 29^ gl i,9''-9ym + 2{fig^-9g2f9^ 

2 2m mii-9^ 2(m// - ^2)2 



-2P1C1 ^1 



m6'^2 - 6*2^1 
m/j, — ^2 



PiM ^2 



m9gi - 6*2^2 
m/i — ^2 



-2C7iM f "^^^1^2 - /^6>ff^ ^ n9{mg2 - 9gi 

\ m/j, — 6*2 2{mn — 6*2)2 j-,^^ _ ^2^2 



m 



2 ^ 1 ^ ii9^ - 9'^/m + 2{iigi-9g2f9^ 



2m 



-m9 



mngig2 - 1^9 gf 



9 



2(m// - 6*2)2 
1^9 {mg2 -9giY 



mfJL — ^2 



2 



2{mn - 6*2) (m// - ^2)2 

32^ ^/l2 M 



^2 gi 



m9g2 - 6*2^1 
m/x — ^2 

2/32- 



2 1 2m%i^2-2g2^i mg2_g4/^^2(mgr2-%i)"g 
2// m// - ^2 + 2(m/x - ^2)2 
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2 V 2// 



+ 



-2P2C.(^.- ^^^^-^> )-P.M(ff,- 



-2C2M 



m// — ^2 
/ ixmgig2 - n9gl _ 



2{mii - ^2)2 
m9g2 - d^gi 



1^0 {mg2 -9giY 

2{miJ-9'^) {m^i-9^Y 



9 



2^ 1 _2m9gig2-29'^gl ^m9'^ -9^/^ + @{mg2-9gif9'^ 



mil — 6*2 



6* 



2(m// - ^2)2 
(4^f + 6^i/;u)(m^2 -^^1) 



m// — 6*2 



36*^ 



[m^2 - , 3 [m^2 - Ogi\ 



+ 



[m// - ^2] 3 ^ 2// [m// - ^2] 2 



+ 



I g4 / ["•^^2 - OgiY' ^ 3 [mff2 - Ogif ^ 3 



[m// - ^2]4 ^ _ ^2]3 4^2 _ ^2]2 



m 



^, 2 1 2^%ig2 - 29''gl ii9^ - 9^/m + 2(/xgi - ^^2)'^ ' ' 
2 2m m/x - ^2 + 2(m/x - ^2)2 



— 6£ 



1 

2/i 



li^{mg2 - 9g] 



29 



2(m/i-6'2) (m^u- 6*2)2 

3r^^ Z)^ 13 Q ,.2 



_ lJ^^[mg2 - 9g{f' 3 n'^[mg2 - 9gi 
//^'V [m/i-6'2]3 +2 [m//-^2]2 



+ 



6*2 ^ //[m5f2 -9gi]'^ _^ ^IJ'^[mg2 - 9gi 



|2 3 

+ T- 



[m// — ^2]3 4[m// — ^2]; 



where ^ = {9'^[iig\ + m^r^] ~ 2m/i^i^2)/("7.At - ^^). 



(SI) 
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